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Abstract. We analyze the coupled non-linear integro-differential equations whose solution is the thermody- 
namical limit of the empirical correlation and response functions in the Langevin dynamics for spherical p-spin 
disordered mean-field models. We provide a mathematically rigorous derivation of their FDT solution (for the 
high temperature regime) and of certain key properties of this solution, which are in agreement with earlier 
derivations based on physical grounds. 



1. Introduction 

The complex long time behavior predicted for the thermodynamical limits of a wide class of Markovian 
dynamics with random interactions, is among the fascinating aspects of out of equilibrium statistical physics 
(for a good survey on phenomena such as aging, memory, rejuvenation, and violation of the Fluctuation- 
Dissipation Theorem (FDT), see [7, 8]). This work is concerned with the long time behavior of a complex 
system composed of N Langevin particles xt = (xl)i<i<N € K , each evolving in E and interacting with the 
others through a random potential. More precisely, one considers a diffusion of the form 

(1.1) dx t =-f'(\\x t \\ 2 /N)x t dt- pVHj(x t )dt + dB t , 

where B t is a A^-dimensional Brownian motion, ||x|| denotes the Euclidean norm of x £ M. N and / is a convex 
function. Such models are called spherical (for well chosen /, xt is restricted to stay on a sphere). The mixed 
p-spin, p < m, potential Hj : M. N — ► R is given by 

m 

(1.2) Hj(x)=J2^J E W 1 - 1 *'. 

P= 2 l<ii<— <i p <N 

where the coupling constants J^.-.i are assumed to be independent centered Gaussian variables. The variance 
of Jn...i p is c({ii, . . . , i p })N~ p+1 , where 



tfe! • 



(1.3) c({h,...,i p }) = ]Jk 

k 

and (li, I2, ■ ■ ■) are the multiplicities of the different elements of the set {i\, . . . , i p } (for example, c = 1 when 
ij 7^ iji for any j 7^ j' , while c — pi when all ij values are the same). 

When m = 2, one gets the so-called Sherrington- Kirckpatrick spherical spin glass, which has been studied in 
details in [3]. Given a realization of the coupling constants, the dynamics of (1.1) is invariant for the (random) 
Gibbs measure 

N 

(1.4) fi J N (dx) = ZJ^eM-NfiWx^/N) ~2(3Hj(x))'[[dx\ 



i=i 
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Similar random measures have been extensively studied in mathematics and physics during the last two decades 
(see e.g. [20], for the rigorous analysis of the asymptotics of the free energy of the measure with a hard spherical 
constraint, corresponding to spins on the sphere ||x|| 2 = N). Here, we shall be concerned with the statistical 
properties of the dynamics at high temperature. The natural quantity of interest is the empirical covariance 
function 

1 N 

(1.5) C N (s,t) = -J2 X >1' s ^t> 

i=i 

in the large N limit, and for large t and s. It turns out that the asymptotic behavior of (1.5) strongly depends 
on the way t and s tend to infinity, at least at low temperature. This is a trace of aging: the older it gets, 
the longer the system will take to forget its age. This innocent looking notion of aging is related to deep 
mathematical problems (see e.g. the survey in [13]), and leads to interesting mathematical scenarios, like 
the ultrametric property of the covariance function at low temperature (see e.g. [10]). In [3], the authors 
present a detailed analysis of the aging properties of (1.5) in the special case m = 2, using integro differential 
equations involving the almost sure limit C(s,t) = limjv^oo Cn{s, t). When m ^ 2, closed equations for C 
are obtained in [10] (for the hard spherical constraint) and rigorously derived in [5], where they are called 
Cugliandolo-Kurchan equations. These equations involve also the limit of the integrated response function 

N 



(1-6) XN(s,t) = ^Y,x\Bl 



N 

i=l 

According to [5], fixing T < oo, the random functions Cn and xn converge uniformly on [0, T] 2 , almost surely 
and in L x to non-random functions C(s, t) and x(s, t) = J Q R(s, u)du with R(s, t) — when t > s, R(s, s) = 1, 
and, for s > t, the absolutely continuous functions C, R and K(s) = C(s, s) are the unique solutions in the 
space of bounded, continuous functions, of the non linear integro-differential equations 

(1.7) d s R(s, t) = -f(K(s))R(s, t)+[3 2 J R ( u > t) R ( s > u)v"{C{s, u))du, 

(1.8) d s C{s,t) = -f(K(s))C(s,t) + f3 2 f C{u,t)R{s,u)v"{C{s,u))du + (3 2 f u'(C{s,u))R{t,u)du, 

Jo Jo 

(1.9) d s K(s) = -2f'(K(s))K(s) + 1 + 2[3 2 f ip(C(s, u))R(s, u)du, 

Jo 

where ip(r) = v'(r) + rv"(r), 

m a 2 

(1.10) I/ ( r )=V^r p , m>2, a m £ 0, 

and the initial condition K(0) = C(0,0) > is given. It was shown in Theorem 1.2 of [5] that for 

(1.11) f(r) := f L ( r )=L(r-l) 2 + ^-r 2k , k > m/4, k e L > , 

4fc 

these equations admit a unique solution C(s, t) = C(i, s), R(s, t) and K(s) = C(s, s) in the space of absolutely 
continuous functions on {(s, t) : < t < s}. 1 

Whereas [5] rigorously derives these equations using stochastic calculus and concentration inequalities, 
in the physics literature they are attained via the so-called Martin-Siggia-Rose formalism (see e.g. [6] or 
[9]). An alternative to the latter is to expand the stochastic process Xt perturbatively using diagrams, as 
explained for example in [6] and [8], then average over the disorder the product x s Xt to get equations relating 
the covariance and response functions. The generic form of the family of diagrams one uses in this process 



1 More general choice of diffcrcntiablc /(■) is allowed in [5], but we focus here on the collection given by (1.11). 
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indicates that (1.7)-(1.9) are equivalent to a class of mode coupling equations (see e.g. [6], [8] or [17]). Mode 
coupling approximations were developed in physics to study nonlinear random dynamical systems occurring 
in many contexts like plasma physics, kinetic theory of classical liquids or glasses (see e.g. [15] and the 
references therein) . This method considers a perturbative expansion as a series indexed by diagrams containing 
information on the nonlincarity. The series is then renormalizcd, to produce self-consistent moment equations, 
called mode coupling equations. Kraichnan [16] [17] developed such approximations (called direct interaction 
in fluid mechanics) when considering the solutions to the Navier-Stokes equations in Fourier space under 
random initial conditions, and gave a procedure to perform the statistical closure of the moment equations: 
usually, there is a cascade of moments meaning that the time derivative of the second moment involve the 
third moment and so on. One then looks for good approximations leading to self-consistent equations (c.f. the 
texts [18] and [19], for a variety of statistical closure problems and diagrammatic methods in fluid mechanics). 
More recently, these approximations were applied in the study of super-cooled and strongly interacting liquids 
(glasses), producing very accurate quantitative predictions (see [14], [15]). The fundamental object of interest 
in structural glasses is a correlation function n(t), for which the mode coupling equation is 



where k{r) = Y^™=i \ rP j f° r some constants f2, v and {b p ,p — 1, . . . , to} (see e.g. [14]). The intcgro-differential 
equation (1.12) is similar to what one gets when postulating the FDT ansatz, whereby the solution (R,C) of 
(1.7)-(1.9) is translation invariant, with K(s) constant and R proportional to the derivative of C (for example, 
see (5.9) in the sequel). Of course, the first task when dealing with (1.7)-(1.9) is to show the validity of this 
ansatz, at least for j3 small enough, as we do in Theorem 1.3. Further, the genericity of (1.12) implies that 
any relevant information about its solution is of interest (and in this context see Proposition 1.4). 

The asymptotic behavior of C(s, t) and R(s, t) for large values of t and s is difficult to pin down; in [10], the 
authors propose various scenarios, but no complete description of these asymptotics could be given (see also 
[7] or [8]). A first regime of interest is the so-called FDT regime, in which the fluctuation dissipation theorem 
of statistical physics is expected to hold. In this regime, the covariance should be stationary, that is, for fixed 
s — t = t and t large, the covariance C(s, t) should be approximated well by some function CfdtOr). Further, 
from [12] we know that the response function R(s,t) is then well approximated by some function -RfdtM- In 
this regime we further expect the FDT relation i?fdtC?~) = — 2C f ' dt (r) to hold (the proportionality constant has 
to do with the scaling we employed in the definition of R(s,t) and not with its physical meaning). We shall 
see in this work that the FDT regime holds for small enough /? in which case both Cfdt and i?fdt decay to zero 
exponentially fast. The aging regime is expected to be characterized by covariances and response functions of 
the generic form C(s,t) = C ag i ng (h(t)/h(s)) and R(s,t) = h'(t)h(s)~ 1 R ag in g {h{t)/h{s)) with sub-exponential 
growth of the monotone function h(-) and a polynomial decay to zero and to appear only for (3 > (3 C , the 
dynamical phase transition point of (1.1). These different scenarios are examined in [12], where fixing the 
asymptotic behavior of C according to the above choices Cfdt or C ag i ng , the authors study the solution R of 
the equation 



(1.12) 





where k(s,u) stands either for f} 2 v" '(Cfdt(s — u)) or for [3 2 v" ' (C ag i ng (h(u) /h(s)). Let 




Then it is easy to check that H solves the equation 



(1.13) 
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The main technique of [12] is a formal solution to (1.13) involving non-crossing involutions. We shall also use 
this formula in the sequel to prove various statements concerning the solution C and R of the complete system 
of equations (for example, their exponential decay at high enough temperature) . 

We first show that in the limit L — ► oo, the solutions of the equations (1.7)— (1.9) for our "soft" spherical 
constraint by the potential Jl{ ) coincide with the limiting equations of [10] for the hard spherical constraint 
(where C(2-, 2-) and R(2-, 2-) are considered for v{r) = r p /&). 

Proposition 1.1. For any T < oo, the solution (Rl, Kl,Cl) of the system (1.7)-(1.9) with potential /l(-) 
as in (1-11) and initial condition Kl(0) = 1, converges as L — > oo, uniformly in [0, T] 2 , towards the triplet 
(R, K, C) such that C(t, t) = K(t) = 1 for all t > 0, R(s, t) = for all s <t, and for all s>t, 

(1.14) d s R{s, t) = -n(s)R(s, t)+/3 2 J R(u, t)R(s, u)v"{C{s, u))du, 

(1.15) d s C{s,t) = -n(s)C(s,t) + /3 2 [ C(u,t)R{s,u)v"(C{s,u))du + 1 [ v'{C(s,u))R(t,u)du, 

Jo Jo 

where 

(1.16) n(s) = i + 2f3 2 J ^(C{s, u))R(s, u)du 

Moreover, C(s,t) = C(t,s) is a non-negative definite kernel, with values in [0,1] and R(s,t) > is such that 
R(t, t) = 1 and 

t>t2 

(1.17) 1/ R(s,u)du\ 2 < t 2 - h , 0<ti<t 2 <s<oo. 

Jti 

We prove in Section 3 that for (3 > sufficiently small, the solution (R, C) of (1.14)— (1.16) decays exponen- 
tially fast in s — t. 

Proposition 1.2. There exists f3 > such that for all (3 < [3 there exists Mp < oo and r\p > for which 
the solution of (1.14)-(1.16) satisfies for all s > t > 0, 

(1.18) R(s,t) < e -w(s-t) 

(1.19) C(s,t) < M (i e-^ s - l) . 

Equipped with Proposition 1.2 we prove in Section 4 that for some (3\ > and each j3 < /3±, when s — t = t 
is fixed and t — > oo the solution (R(s , t) , C \s , t)) converges to a limiting pair (i?fdt(r), Cfdt(T)) that satisfies 
the FDT relation iifdt(T) = — 2Cf dt (r) for all t > 0. Upon analyzing in Section 5 the corresponding FDT 
equations (4.15)-(4.17), we establish our main result. 



Theorem 1.3. If < fa. < 1/(2^" (0)) then for any r > 0, 

lim C(r + t,t) =Cfdt(r), 

t— *oo 

and 

lim J2(t + *,*) = flfdt(r) = -2C^ t (r), 

i — >oo 

w/iere Cfdt(-) is the unique [0,l]-valued, continuously differ entiable solution of the equation 

(1.20) £)'(«) = -/ 0(D(u))D'(s-u)dw-6, D(0) = 1 , 

Jo 

/or o = 1/2 and 0(x) = b + 2j3 2 i>'{x). Moreover, both Cfdt(-) and Rmt{-) decay exponentially to zero at infinity. 

A key ingredient of the proof of Theorem 1.3 is the analysis of the equation (1.20), which is of some 
independent interest. Specifically, in Section 5 we prove the following. 
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Proposition 1.4. Suppose b > and <fi : [0, 1] — > R is non- decreasing, of Lipschitz continuous derivative such 
that 

(1.21) sup {>(x)(l ~x)} > b. 

0<x<l 

Then, (1.20) has a unique solution in the space of [0,l]-valued continuously differ entiable functions. This 
solution is twice continuously differ entiable, strictly decreasing and converges for s — > oo to 

Doo := sup{x e [0, 1] : <j>(x)(l - x) > b} . 

Further, in case cp(-) is convex and 0(1) > 2y/b(j>'(l), the derivative of the solution of (1.20) decays exponen- 
tially to zero with some positive exponent. 

Setting b = 0(D OO )(1 — Z?oo) and considering x J. in (1.21), it is not hard to verify that necessarily 
(('{Doc) > 0'(£>oo)(l — -Doc)- In the remark following the proof of Proposition 1.4 we observe that the condition 

(1.22) 0(D OO )>0'(D OO )(1- J D OO ), 

is necessary for the exponential convergence of D'(s) to zero as s — > oo when </>(•) is convex. As (1.22) is easily 
seen to be equivalent to (f)(1) > 2^/&0'(l) when (/>(a;) is linear, we speculate that it actually characterizes the 
exponential convergence of the solution of (1.20). 

Of course, in case of <(>(•) = b + 20 2 v'(-) and b = 1/2 the condition (1.21) holds (try x — 0). Setting 
C € (0, oo) via 

(1.23) = sup{i/(x)(l - x)x- 1 : x e (0, 1]} , 

it is easy to check that in this case Doc = if < C whereas Z?oc > for > C . Further, considering x — > 
in (1.23) wc find that 1/(4/3;?) > i/'(Q), so the condition (1.22) then holds for any < C . This indicates that 
though the values of 0o and 0\ in our proofs of Proposition 1.2 and Theorem 1.3 arc quite small, both should 
match the predicted dynamical phase transition point C of our model. 

Indeed, subject to a heuristic ansatz, similar to what can be found in the physics literature, we show in 
Section 6 why the equation (1.20) with 6=1/2 and <f>(x) = 7 + 20 2 v'(x) describes the FDT regime at all 
temperatures, albeit in general with 7 = ^(0) ^ b (see Proposition 6.1). More precisely, the physics prediction 
for the choice of 7 is as follows. Let x* denote the largest value of x < 1 attaining the supremum in (1.23) 
of the polynomial h(x) = v'(x)(l — x)/x (with h(0) = v"(0)). Setting g(x) = v"(x)(l — x) 2 we have that 
g(x*) = h(x*) = l/(A0 2 ) as a consequence of the optimality condition h'(x*) = if x* > (while trivially 
g(0) = h(0)). From the continuity of g(-) and the fact that g(l) = it follows that 

(1.24) q{0) := sup{x € [0, 1] : A0 2 v"{x){l - x) 2 > 1} , 

is in [a;*, 1) for all > C , with q(0) > x* as soon as > C . Setting ~/{0) = 20 2 [v" [x){l - x) - v'(x)] for 
x = q(0) we find that the condition (1.21) of Proposition 1.4 applies also for > C . Further, the expected 
limit Dqo of Cfdt(T) as r — > 00 is the preceding q(0) which is strictly positive for > C , indicating the onset of 
the aging regime at C . Another indication of the onset of aging is the fact that when > C the exponential 
convergence of Cfdt(-) and -Rfdt(-) is lost (i.e. the preceding choice of j(0) leads to equality in (1.22)). The 
physics prediction suggests also that 

I 7 (/3) = 2 / 1p(C 'aging (A)) Raging (A)rfA > , 

Jo 

where I 7 (/3) = l{0) — \ + 20 2 q(0) l/ ' {q{0)) is strictly positive when q(0) > 0, and in particular, whenever 
> 0c (c.f. the remark following Proposition 6.1). The precise nature of the dynamical phase transition at 
0c depends on whether q(0 c ) is strictly positive (as is the case for example whenever = 0), or not. The 
physics ansatz of one aging regime with Raging (A) = AC' aging (\) and C ag ing{') monotone increasing on [0, 1] 
from C a gi„g (0) = to C ag i ng (l) = q(0), proposed in [10] for the case in which q(0 c ) > 0, thus allows one to 
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set the positive constant A = I 1 rm/(f3 2 q(f3)u'(q(f3)) in terms of the FDT solution. Unfortunately, even under 
this ansatz, the form of C ag i ng {-) is yet unclear. 

Finally, note that it is crucial for our analysis to have v'(0) = 0, an assumption we make throughout this 
work. That is, our analysis holds in the absence of a random magnetic field. 

2. Limiting exactly spherical dynamics 

Proposition 1.3 and (2.13) of [5] show that sup t>0 K(t) < oo. Further, as C(s,t) is the limit of the 
empirical correlation functions Cj\r(s,t) it is a non-negative definite kernel on M + x K + and in particular, 
C(s, t) 2 < K(s)K(t), whereas since x(s, t) = J * R(s, u)du is the limit of Xn(s, t), it follows from the definition 
of xn(s, t) that 

(2.1) | / R(s, u)du\ 2 < K(s)(t 2 - ii) , < h < t 2 < s < oo . 

To complete the proof of Proposition 1.1, we first prove that any solution (R, C, K) of (1.7)— (1.9) consists 
of positive functions, a key fact in our forthcoming analysis. 

Lemma 2.1. For any f : M + — ► K whose derivative is bounded above on compact intervals and any K(0) > 0, 
a solution (R,C,K) to (1.7)-(1.9), if it exists, is positive at all times. 

Proof: By definition K{t) > for all t G K + . Suppose that 

S = inf{u > : C(t, u) < for some t < u} < oo . 

By continuity of (C,K), since K(0) > also S > 0. Set A(s,t) = cxp(- n{u)du) > for fj,(u) = f'(K(u)) 
which is bounded above on compact intervals, and R(s,t) — A(s,t)H(s,t). Then, by [12], for s > t, 

. 2n 

(2.2) ff( S ,t)=i+^/3 2 " yi / n ""mut^Udtj 

™>! <7eNC„ -* 1 "- t2n - s *ecr(<r) 3=1 

where NC„ denotes the set of involutions of {1, • • • , 2n} without fixed points and without crossings and cr(er) 
is defined to be the set of indices 1 < i < 2n such that i < a(i). Consequently, 

R(s, t) > A(s, t) > for t < s < S , 

and thus, (1.8) implies that also 

C(s, t) > K(t)A(s, t)fovt<s<S. 

Note that in the last two estimates we used the fact that v'(-) and v"{-) are non negative on R + . Similarly, 
from the equation (1.9) we see that d s [A(s, 0)~ 2 K (s)] > A(s,0) -2 for all s < S resulting with 

K(s) >K{0)A(s,0) 2 + f A(s,v) 2 dv>0 
Jo 

Hence, the continuous functions R(s,t),C(s,t) and K{s) are bounded below by a strictly positive constant 
for t < s < S in contradiction with the definition of S. We thus deduce that 5* = oo and by the preceding 
argument both R(s,t) and C(s,t) are positive functions. □ 
We next show that if (Rl,Cl, Kl) are solutions of the system (1.7)— (1.9) with potential /l(-) as in (1.11), 
then Kl(s) — > 1 as L — > oo, uniformly over compact intervals. Specifically, 

Lemma 2.2. Assuming Kl(0) = 1, we have that Kl(s) > 1 for all L > and s > 0. Further, for any T 
finite there exists B(T) < oo, such that Kl(s) < 1 + _B(T)L _1 for all s <T and L > B(T). 
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Proof: We first deal with the lower bound on Kl(-). To this end, fix L > and let g(x) := 1 — 2xf' L (x). It 
is easy to check that g(x) = 1 + ALx(l — x) — x 2k > for x e [0, 1] and g(x) < for x > 1. By Lemma 2.1, 
we know that the functions Rl(-, ■)> an d Cl(-, ■) are non negative, as is ^>(r) for r > 0, so from (1.9) we get 
the lower bound d s Ki J (s) > g(Ki J (s)). Thus, with <fi(x) a differentiable function that is strictly increasing on 
[0. 1] and such that <j>(x) = 1 for all x > 1, we find that 

d s ^(K L (s)) > cf>'(K L (s))g(K L (s)) > 0. 

Consequently, <P(Kl{s)) > 4>{Kl{Q)) — <p(l) = 1 for all s > 0, implying by the choice of </>(•) that Kl(s) > 1 
for all s > 0. 

Turning now to the complementary upper bound, recall that ip(r) is a polynomial of degree m — 1, hence 
there exists k < oo such that ip(ab) < k(1 + a 2 )" l / 2 (l + 6 2 ) m / 2 for all a, b. Thus, by (2.1), the monotonicity of 
ij)(r) on R + and the non-negative definiteness of Cl(s, u) we have that for any s,t,u> 0, 

^(C l {s,u))<k(1 + K l {u))^{1 + K l (s))^, f R L {s,u)du<^tK L {s), 

Jo 

and from (1.9) we find that 

(2.3) d s K L {s) < g{K L {s)) + 2/3 2 k(1 + su P K L (u)) m y/K L (sj^ . 

U<S 

Setting now B(T) = 1 + 4(3 2 K3 m VT and fixing T < oo and L > B(T), let 

r := inf{u > : K L (u) > 1 + BL^ 1 } . 

By the continuity of Kl{-) and the fact that Kl(0) = 1 < l + BL^ 1 , we have that r > and further, if r < oo 
then necessarily 

K l {t) = sup K l (u) = 1 + BL- 1 < 2. 

U<T 

Recall that g(x) < 1 + 4L(1 — x) when x > 1, whereas from (2.3) we see that if r < oo then 

9 s X L (s) u=T < l-4S + 4/3 2 K3 m V7. 

Recall the definition of t < oo implying that d s KL(s) > at s = t. Hence, our choice of B = B{T) results 
with t > T. That is, K L (s) < 1 + BL^ 1 for all s < T and L > B(T), as claimed. □ 
Let Hl{s) = f' L (KL(s)) and /il(s) = d s KL(s). Fixing hereafter T < oo, we next prove the equicontinuity 
and uniform boundedness of (Cl, Rl, Kl 1 Ml, ^l), en-route to having limit points for (Cl, Rl,Kl)- 

Lemma 2.3. The continuous functions (Cl(s, t), Rl(s, t)) and their derivatives are bounded uniformly in 
L > B(T) and < t < s < T. Further, the continuous functions (^l(s), /il(s)) and their derivatives are 
bounded uniformly in L > B(T) and s G [0, T]. 

Proof: Recall that by Lemma 2.2, for any L > B(T), 

(2.4) sup\K L (s)-l\<^p- . 

s<T L 

Consequently, the collection {CL(s,t),0 < t < s < T, L > B} is uniformly bounded. By (2.4) and our choice 
of fh if) ! we have that 

\li L (a)\ < 2L\K L (s) - 1| + K^s)^- 1 < 2B(T) + 2 2k ~ 1 , VL > B(T), s < T , 

whereas by (2.2) the collection {R L (s, t), < t < s < T, L > B} is also uniformly bounded. Further, since 

(2.5) h L {s) = l-2K L {s)ii L {s) + 2(3 2 f ^(C L (s,u))R L (s,u)du , 
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it follows from the uniform boundedness of Kl, \il, Cl and Rl that {/il(s), s £ [0, T], L > B} is also uniformly 
bounded. By the same reasoning, from (1.7) and (1.8) we deduce that d s CL(s,t) and d s B,L(s,t) are bounded 
uniformly in L > B and s,t £ [0, T]. 

Next, differentiating the identity (2.2) with respect to t, we get for / = /j, that 

„ In 

d t H L ( S ,t)=Y,p n e / n ""(cutut^Udtj, 



where NC n denotes the finite set of non-crossing involutions of {1, . . . , 2n} without fixed points. With the 
Catalan number |iVC„| bounded by 4", and since Cl^, € [0,2] for ti,t a ^ < T, L > B(T), we thus 
deduce by the monotonicity of r i— > v" (r) that 

< d t H L (s, i) < £ 7^T)! 4V "( 2 ) n ( s - < ) 2 "" 1 - 

n>l 

so d t H L (s,t) is finite and bounded uniformly in L > B(T) and < t < s < T. Since 

d t R L (s,t) - Li L (t)R L (s,t) - e--l t s ^ du d t H L (s,t) , 

we thus have that \d t RL(s, t)\ is also bounded uniformly in L > B(T) and < t < s < T. 

In the course of proving [5, Lemma 4.1], one finds that d t C(s, t) = R(s, t) + D(s, t) for the function D(s, t) 
of [5, (4.2)]. Consequently, in case of / = /l we have that 

d t C L (s,t) = R L (s,t) - n L (t)C L (s,t) + f3 2 f C L (s,u)R L (t, u y'(C L (t,u))du + p 2 f v' (C L (t,u))R L (s,u)du , 

Jo Jo 

which of course is also bounded uniformly in L > B(T) and < t < s < T. 

Turning to deal with /il(-), setting gL{r) ■= [,/^(r)r]' — 2L = AL{r — 1) + fcr 2fe_1 , we deduce from (2.4) 

that \g L {K L {s))\ < AB{T) + k2 2k for any s < T and L > B(T). Differentiating (2.5) we find that d s h L (s) = 

-ALh L {s) + k l (s) for 

k l {s) = -2g L {K L {s))h L {s) + 2f3 2 d s f ^{C L {s,u))R L {s,u)du , 

Jo 

which is thus bounded uniformly in L > B(T) and s < T (in view of the uniform boundedness of h,L, Cl, 
Rl, 8 s Cl and 8 s Rl)- Further, recall that Kl(0) = 1, so by (1.9) and our choice of /i,(r) we have that 
h L (0) = 1 - 2f' L (l) = 0, resulting with 



h L (s) = f e - iL{s ~ u) K L (u)du . 
Jo 



Since A(T) = sup{|k l (m)| : L > B(T),u< T} is finite, we thus deduce that for L > B(T), 
(2.6) sup \h L (a)\ <4P' 

S<T 4iy 

and the uniform boundedness of |9 s /il(s)| follows. 

Finally, by definition, 9 s /x_l(s) = f'l(KL(s))hL(s), yielding for our choice of fa that 

\d,n L (a)\ < (2L + (2k - l)2 2k - 3 )\h L (s)\ , VL > B(T),s < T , 

which by (2.6) provides the uniform boundedness of |9 s /xl(s)|. □ 

Proof of Proposition 1.1. By Lemma 2.3 we have that (Cl(s, t), Rl{s, t)), L > B(T) are equicontinuous 
and uniformly bounded on < t < s < T. Further, [K Ll n L ,h L ) are then equicontinuous and uniformly 
bounded on [0, T\. By the Arzela-Ascoli theorem, the collection (Cl, Rl, Kl, I^l, ^l) thus has a limit point 
(C, R, K, fi, h) with respect to uniform convergence on [0, T] (or < t < s < T, whichever is relevant). 
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By Lemma 2.2 we know that the limit K(s) = 1 for all s < T, whereas by (2.6) we have that h(s) = for 
all s <T. Consequently, considering (2.5) for the subsequence L n — > oo for which (Cx n , i?L„, ^Cl„j Mini ^l„) 
converges to (C, i?, X, /U, h) we find that the latter must satisfy (1.16). Further, since iiz,(t, i) = 1 and Cx(i, t) = 
X L (t), integrating (1.7) and (1.8) we find that R L (s, t) = 1+]*/ A L {6, t)dO and C L (s, t) = K L (t)+ f° B L (0, t)d6, 
for 

A L (9,t) = -ii L (e)R L {6,t)+p 2 J R L (u,t)R L (6,uy'(C L (6,u))du, 

B L (9,t) = -n L (8)C L (9,t)+/3 2 [ C L (u,t)R L (9,uy'(C L (9,u))du + (3 2 [ v> \C L {9,u))R L {t,u)du . 



Since A^ n (s, t) and Bl„ (s, t) converge uniformly onO<t<s<Tto the right-hand-sides of (1.14) and (1.15), 
respectively, we deduce that for each limit point (C, i?, /i), the functions C(s, t) and R(s, t) arc diffcrentiable in 
sonO<t<s<T and all limit points satisfy the equations (1.14)— (1.16). Further, Cl(s, t) are non-negative 
functions, and also symmetric non-negative definite kernels with Cl(4, t) — > 1. Consequently, each of their limit 
points corresponds to a [0, 1] -valued symmetric non-negative kernel on [0, T} 2 . Similarly, since R^t.t) = 1 
and Rl(s, t) satisfy (2.1), the same applies for any limit point R(s, t). The latter are extended to functions on 
[0, T} 2 by setting R(s, t) = R L (s, t) = whenever s < t. 

Finally, it is not hard to verify that the system of equations (1.14)— (1.16) with C(s,t) — C(t,s), C(t,t) — 
R(t, t) = 1 and R(s, t) — for s < t, admits at most one bounded solution (R, C) on [0, T] 2 . Indeed, considering 
the difference between the integrated form of (1.14)— (1.15) for two such solutions (C,R) and (C,R), since v" 
is uniformly Lipschitz on [0, 1], the functions AR(s,t) = \R(s, t) — R{s,t)\ and AC(s,t) — \C(s,t) — C(s,t)\ — 
AC(t, s) are such that for < t < s < T, 

(2.7) AR(s,t) < Ki[f AR(v,t)dv+ [ h(v)dv] 

(2.8) AC(s,t) < AC(v,t)dv + h(t)+ ( h{v)dv] 

Jt Jt 

where h(v) := J^[AR(v,u) + AC(v,u)]du and K\ < oo depends on T, /?, v(-) and the maximum of \R\, 
\C\, \R\ and \C\ on [0,T] 2 . Integrating these inequalities over t G [0, s], since AR(v,u) = for u > v and 
AC(v,u) — AC(u,v), we find that 



AR(s,t)dt < n 2 / h(v)dv, 
Jo 

/ AC(s,t)dt < n 2 / h(v)dv, 
Jo Jo 

for some finite constant k 2 (of the same type of dependence as ki). Summing the latter two inequalities we 
see that for all s€ [0, T], 



< h(s) < 2k 2 [ h(v)dv. 
Jo 



Further, h(0) = 0, so by Gronwall's lemma h(s) = for all s e [0,T]. Plugging this result back into (2.7)-(2.8) 
and observing that AR(t,t) = AC(t,t) = 0, we deduce that AR(s,t) = AC(s,t) = for all < t < s < T, 
yielding the stated uniqueness. 

In conclusion, when L — > oo the collection (Cl, Rl, Kl) converges towards the unique solution (C,R,K) 
of (1.14)-(1.16), as claimed. □ 

3. Exponential decay for small values of (3 

We consider here the solution (R, C) of (1.14)-(1.16) and prove Proposition 1.2 about its exponential decay 
in s — t for all j3 > sufficiently small. 
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3.1. Exponential decay of R(s,t) for (i small. 



Lemma 3.1. If (3 < l/(4yV'(l)) then R(s,t) < e^ 8 "*) for 8p = \- 2f3y/v"{\) > and all s>t. 

Proof: By [12], for s > t, the solution (R,C) of (1.14)-(1.15) is such that R(s,t) = A(s,t)H(s,t) for A(s,t) = 
exp(— J t s (j,(u)du) with /j,(u) of (1.16) and H(s,t) of (2.2). Further, since C(s,t) £ [0,1] for all s,t > 0, it 
follows that v"(C(s,t)) < < oo. Recall that |jVC n | < 4", hence we deduce that 

H(s, t) < 1 + Y, /3 2 "4V"(1)" (S ^ )2 " < ^(2/3v^(l)) fc ^i^ = e 2/3\^ 77 W(-*) . 

n>l ' fe>0 

By Lemma 2.1 we know that C and R are non-negative functions, hence fi(u) > \ (by (1.16)), resulting with 
A(s,t) < e~ (s ~*)/ 2 . In conclusion, 



R(s,t) = A(s,t)H(s,t) < e 



(2/? x / I /"(l)-i)( 5 -t) = -5 3 (s-t) 



where Sp > for (3 < l/(4yV 7 (T)), as claimed. □ 
3.2. Exponential decay of C(s,t) for /3 small. 

Lemma 3.2. For some /3 > and any /3 < [3q there exist Mp < oo and < np < 5p such that 

C(s,t) < Mpe-^ 8 ^ . 
Proof: From equation (1.8) we get that for s > t > 0, 

(3.1) C(s,t) = A(s,t) + f3 2 ( A(s,v)h(v,t)dv + f3 2 [ A(s,v)I 2 (v,t)dv 

Jt Jt 

with A(s, v) = cxp(— f* fi(u)du), 

(3.2) h(v,t) = f C{u,t)R(v,u)v"{C(v,u))du, 

Jo 

(3.3) h{v,t) = I v'{C(v,u))R(t,u)du. 

Jo 

For [3 < l/(4 x /i/"(l)) we know from Lemma 3.1 that R(s,t) < e -^( s -*). With C(v,u) = C(u,v) £ [0,1] and 
since v' (r) < v"(l)r and v"{r) < v"(l) for < r < 1, we get that for w > i, 

Ji(«,t) < / C{u,t)e- 5 ^ v - u Uu, 

Jo 

h(v,t) < ^"(l)(5^ 1 supC(u,w). 

u<t 

Recall that fi(u) > 1/2, so A(s, w) < e~^ s ~ v ^ 2 . Hence, with the symmetric function A(i, s) := sup u<t „ <s C(u, v) 
we deduce from (3.1) that for s > t > 0, 

A(t, s) < e -^ s -*) + /?V'(1) J* e-^ s - v ^ [jf C(u, t)e-^ {v - u) du + 6/ A(t, v)]dv 

< e -K s -*) + /3V'(1) f e ~i( s - v *> f e'^-^dudv 
Jt Jo 

+/3V'(1) jf A(t,^)[^ 1 e-5( S -")+ J e-^ S - V) - 5 ' ,(v - U) du]dv 

It is straightforward to see that for any 5 £ (0, 1/2) and s > t, 



(3.4) 



J e -i(«-«)-«(f-t) dw < 2(1 - 2S)- 1 e - 5(s - t) 
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and with Sp G (0, 1/2) wc thus obtain for s > t the bound 

A(t,s) < Mpe-W"-Q + Ap J A(t,v)e~ s ^ s - v) dv, 

with Mp = 1 + 2/?V'(l)(l - 25p)- 1 8- 1 and Ap = 2/3V'(l)<^ 1 . Therefore, fixing t > 0, the function 
h t (s) = e 5 "( s -*) A(t, s) satisfies 

h t (s) < Mp + A p h t (v)dv, s>t, 

and so by Gronwall's lemma h t (s) < Mpe Af) ( s ~ t \ We therefore conclude that for any s > t, 

C(s,t) < Mpe- {s ?- A ^ s - t] , 

which proves the lemma since for (3 — > we have that Sp — > 1/2 while Ap — ► (and so r\p = Sp — Ap > for 
any j3 > small enough). □ 

4. Getting the FDT equations 
With r = {(s,():0<f<s}cl + xM + ,we consider the map * : (R, C) -> (#, C) on 

A + = {(R, C) e C(r) x C(R+ x R+) : i?(t, t) = C(i, t) = 1, #(s, i) > 0, C(s, i) = C(t, s) > 0} , 
such that for s > i, 

(4.1) <9 s .R(s, i) = -fi RtC (s)R(s, t)+0 2 R(u, t)R(s, u)v"(C(s, u))du, 

(4.2) d s C(s,t) = -/j, RtC (s)C{s,t)+/3 2 [ C(u,t)R(s,uy'(C(s,u))du + /3 2 [ u'(C{s,u))R{t,u)du, 

Jo Ja 

where R(t,t) = C(t,t) = 1, C(t,s) = C(s,t) and 

(4.3) hr, c (s) = \+P 2 J q *KC(s, u))R(s, u)du. 

Assuming (R,C) G A + , we have that hr,c( s ) > 1/2 is continuous and further, by [12] there exists a unique 
non-negative solution R(s,t) of (4.1) which is continuous on T (see for example (2.2) for existence, uniqueness 
and non-negativity of the solution, and the proof of Lemma 2.3 for the differentiability, hence continuity of 
R(s,t)). With C > and R > 0, clearly there is also a unique non-negative solution C(s,t) to (4.2) which is 
continuous on T and due to the boundary condition C(t,t) = 1, its symmetric extension to M + x R + remains 
continuous. Thus, $(A + ) C A + . 

We proceed to show that for small /3 and a suitable choice of the positive constants 5, r, p, c the solution 
(R, C) of (1.14)- (1.15) is a fixed point of the mapping * on the space 

S(5,r,p,c) = {(R,C) G A(5,r,p,c) : 3i? fdt (r) = lim R{t+r,t) , 3C7 fdt (-r) = C7 fdt (r) = lim C(t+T,t) , Vr > 0} 
where 

A(S,r,p,c) = {(R,C) G .4+ : C(s,t) < ce' 5 ^, R(s,t) < p(r\s - t\ + i)-3/a e -«(»-t) for all a > t j. 
This of course implies that the solution (R, C) of (1.14)- (1.15) is such that the functions 

(4.4) R m (r) = lim R(t + T,t), 

t— >oo 

(4.5) C fdt (r) = lim C(t + T,t), 

t— >oo 
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are well defined for all r > 0. Further, for any (R,C) G S(S, r, p, c) the corresponding functions (i2fdt,Cfdt) 
are clearly in the set 

B(S, r, p, c) := {{R, C) g B(R+) x : (7(0) = R(0) = 1, < C(r) = C(-r) < ce" 5 ^, 

< i?(r) < ^(rT + l)- 3 / 2 e- 5r } , 

so in particular, (4.5) holds for all r G R and (i?fdt, Cfdt) of (4.4)-(4.5) are of exponential decay. 

To this end, we start by finding constants 6, p, c and r = r{(3) for which S{8 1 r 1 p 1 c) is closed under the 
mapping 

Proposition 4.1. There exist finite, positive universal constants c\ and Mi > 2, such that for c = 2, p = c\, 
S = 1/6 and r = P^Jcj < l/(3Mi), both 

(4.6) *(A{5,r,p,c))CA(5,r,p,c), 
and 

(4.7) *(S(5,r,p,c))CS(5,r,p,c). 

Proof of Proposition 4.1: We first verify that (4.6) holds. To this end, setting R(s,t) = A(s,t)H(s,t) for 

(4.8) A(a,t) = e -/>B.c(«)d« ; 

we have that H{t,t) = 1. Further, from [12] we have that for any (s,t) G T, 



^ 2n 

(4.9) H( S ,t)=i+Y,p 2n e / n ""mut^Udtj. 



Consequently, as |NC„| = (2tt) 1 J_ 2 a; 2 "V4 — £ 2 e?x and C(u, u) G [0,c] for (R,C) G ^4(<5, r, p, c) , we have the 
bound 



(4.10) #(M) < 5>V'(c))» E / II 1 

™>o CT£ NC n - 1 -"-* 2n - s J= 1 

= (27T)- 1 J 2 e^V^i-^^^dx. 
It is well known (see for example [3, (3.8)]) that for some universal constant 1 < ci < oo and all 6, 

f-2 



(^r 1 / e e V4-.T 2 dx < ci(l + |0|)- 3/2 



e 2|fl| 



from which we thus deduce that 



(4.11) tf(s, i) < ci(l + 0y/v»(<ij(* - t))~ 3/2 e^V^K'-*) . 

Further, since (R,C) G ^4 + , we know that pn.c{u) > 1/2 resulting with A(s,i) < e~( s ~*)/ 2 . It then follows 
that for our choice of p = Ci, r = f3yjv"(c) and <5 = l/6<l/2 — 2r, 



R{s,t) < ci(l + /3y^^)(s-t))- 3/2 e- ( 2- 2 ^\/^))( s -*) 
< p(l + r( S -t))- 3 / 2 e- 5 ( s -*). 
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Considering next the function C, recall that C(t,t) = 1 and for (R,C) G A(S, r, p, c) we have that 
v'(C(v,u)) < v"(c)C{v,u) and A(s,v) < e~ {s ~ v)/2 for all v < s and u. Thus, we get from (4.2) that for 
(s, t) e r and (R, C) G A(6, r, p, c), 

pS pV pt 

C(s,t) = A(s,t) + fl 2 A(s,v)dv[ C(u,t)R(v,u)v"(C(v,u))du+ v' '(C(v,u))R(t,u)du\ 
Jt Jo Jo 

< A(s,t)+f3 2 cpv"(c) f e-^ s - v ^e-^ v -^dv[j (r{v - u) + l)~Uu + f (r(t - u) + du] 

Jt Jo Jo 

< [1 + 2Kr- 1 2 cpu"{c)A{s, t){^ - d)- 1 ^ 3 -^ 

where in the last inequality we have used (3.4) for S — 1/6 < 1/2 and 

/•OO 

K := / (0 + l)-id6 = 2. 
Jo 

This shows that (R, C) e A(S, r, p, c) since our choices of c = 2, S = 1/6, p = c\ and r = fly/ v"(c), result with 

2Kcpf3 2 v"(c) 
( 5 -S)r 

once we take Mi = AKc\ and fl sufficiently small for r = r(fl) < l/(3Mi). 

Our next task is to verify that (4.7) holds. That is, assuming that (R, C) G S(5, r, p, c) we are to show that 
the limits (Rfdt,C td t) exist for the solution (R,C) of (4.1)-(4.3). To this end, recall that by (4.8), (4.9), and 
(3.1), for any t > and t > v > 0, 

A(t + T,t + V) = e -(r-v)/2 e -0>j : i o (t+u,t)du 
R(t + T,t) = A(t + T,t)H(t + T,t) 

. 2n 

H(t + r,t) = i + ^/3 2n Yl / II y" {C{t + 6^ + 6^)) J{de j 

«>1 CTG NC„ O - 01 -"'- 92 "- T ^cr( ( T) 3 = 1 

PT pt 

C(t + T,t) = A(t + T,t)+f3 2 A(t + T,t + v)h(t + v,t)dv + i3 2 A(t + r, t + v)I 2 (t + v, t)dv 

Jo Jo 

where by (4.3), (3.2) and (3.3), 

I (t + v,t) = J ip(C(t + v,t + u))R(t + v,t + u)du 

/V 
C(t + u, t)R(t + v,t + u)v"(C{t + v,t + u))du 

,o 

I 2 (t + v,t) = J v'(C(t + v,t + u))R(t,t + u)du. 

Since ip(-), v"{-) and v'(-) are continuous and (R,C) G S(5,r, p,c), as t — > oo the bounded integrands in the 
preceding formulas converge pointwise (per fixed u = v — 9) to the corresponding expression for (i?fdt, Cf<jt)- 
Further, by the exponential tails of (R,C) G A(5,r, p,c) the integrals over [— t, — m] in the formulas for Iq, I\ 
and I2, are uniformly in t bounded by pS~ 1 ip(c)e~ Sm . Thus, applying bounded convergence theorem for the 
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integrals over [— m, v], then taking m — > oo, we deduce that for each fixed v > 0, 

/•CO 

Jo := Hm I {t + v,t)= / i>(C{ dt {6))R{ dt (6)d6, 

/•OO 

(4.12) ?!(«) := lim h(t + v,t)= C idt (v - 9)R idt (9y (C fdt (8))d9, 

i^OO 

/•OO 

(4.13) J 2 («) := lim J 2 (t + «,*) = / */(Cf d t(0))i?fdt(0 - v)d0. 

By the preceding discussion we also know that < Ii(t + v, t) < pi/>(c)(5 _1 for 2 = 0,1,2 and all v, t > 0. Thus, 
with A(t + t, t + v) e [0, 1], by bounded convergence for each r > v > 0, 

A(t-v) := lim A(i + r,i + t>) = e -(^)(i/2+/3 2 ?o) _ 

i — >oo 

C fdt (r) := lim C(i + r, t) = A(r) + [3 2 f A(r - v)h(v)dv + 1 f A(t - v)T 2 (v)dv . 
We also have that for any n € N, all a e NC„ and each fixed 0i, . . . , 6* 2 „ > 0, 

lim II v"{C{t + e i ,t + e a(i) ))= [] ^'(Cm(0i-0a(i))), 

By bounded convergence, we have the convergence of the corresponding integrals over < 0\ < ■ ■ ■ < 9 2n < t. 
Further, the non-negative series (4.9) is dominated in t by a summable series (see (4.10)), so by dominated 
convergence, 

,• 2n 

H {dt (T):=limH(t + T,t) = l + J2P 2n E / II ""(^(ft - ^(»))) H d6j . 

«>1 CTG NC„ °- ei -'"- 02 "- r iecr( ( T) j=l 

It thus follows that 

R ldt (T) := lim R(t + r, t) = A( T )H m (r) , 
exists for each r > 0, which establishes our claim (4.7) (we have already shown that Cf dt (T) exists). □ 

We next show that ^ is a contraction on S(S, r, p, c) and provide the set of equations that characterizes the 
functions Rf dt and Cfdt of Theorem 1.3. 

Proposition 4.2. For S, p, c and r — r{[3) of Proposition 4-1, if is small enough then ^ is a contraction on 
A(S, r, p, c) equipped with the norm 

(4.14) \\(R,C)\\ = sup \R(s,t)\ + sup \C(s,t)\, 

(s,f)er s,teR+ 

and the solution (R,C) of (1.14)-(1.15) is also the unique fixed point of ^ in S(S,r,p,c). Consequently, the 
functions (-Rfdt, Cfdt) of (4.4) -(4.5) are then the unique solution in B(S,r,p,c) of the FDT equations 



(4.15) 


R'(t) 


(4.16) 


C'(r) 


(4.17) 


M 


for I = 


0. 



-pR(r) + (3 2 / R{T-6)R{6)v"(C{6))d6, 
Jo 

p oo /*oo 

-hC{t) + f3 2 C(t - 9)R{8)v"{C{8))d9 + f3 2 v'{C{9))R{8 - r)d9 + I, 

1 

2 



\+p 2 I iP(C(9))R(9)d9 + I, 
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Proof of Proposition 4.2: Keeping S, p, c and r = r(/3) as in Proposition 4.1, we first check that for any 
(3 small enough, ^ is a contraction on A(S, r, p, c) equipped with the uniform norm ||((_R, C)|| of (4.14). To 
this end, we consider the pairs (Ri,Ci) = ^(Ri,Ci) for (Ri,d) £ A(6,r, p,c), i = 1,2. Denoting hereafter in 
short A/(s,t) = |/i(s,i) — /2(s,i)| and A/(s) = sup M< „ <s Af(v,u) when / is one of the functions of interest 
to us, such as R, C, A, H, R or C, we shall show that there exist finite constants Mr — Mr(5,p,c) and 
Mq — Mc{5 1 p, c) such that for any finite s, 

(4.18) AR(s) < M R f3 2 [AR(s) + AC(s)}, 

(4.19) AC(s) < M c p 2 [AR(s) + AC(s)}. 

So, if (3 is small enough for M R f3 2 < 1/3 and M c /3 2 < 1/3, then from (4.18) and (4.19) we deduce that 

- (R 2 ,C 2 )\\ = su P Ai?( s ) + su P AC( s ) < ^[sup AR(s) + sup AC(s)} = - (R 2 ,C 2 )\\ . 

s>0 s>0 o s >0 s >0 J 

In conclusion, the mapping ^ is then a contraction on *4(<5, r, p, c) , since 

(4.20) \\^(R 1 ,C 1 )-^(R 2 ,C 2 )\\ < ^(R^C^-iR^C^W, 

for any (Ri,Ci) G .4.(5, r, p, c), i = 1, 2. 

The challenge in deriving (4.18)^(4.19) is to get bounds that are uniform over (s,t) e T. In doing so, we 
use the tail estimates for (Ri(s,t),Ci(s,t)), i = 1,2 (which hold for all functions in A(S 7 r, p, c)), in order to 
improve upon the arguments of [5, proof of Proposition 4.2], where the Lipschitz bounds are derived for finite 
time intervals. 

• The Lipschitz bound (4.18) on R. 

We rely on the formulas (4.9) and R(s, t) = H(s, t)A(s, t). Indeed, since C\ and C 2 are [0, c]-valued symmetric 
functions, ti e [0, s] and both v"(-) and v'"(-) are non-negative and monotone non-decreasing, it follows that 
for any n, t 2n < s and a <G NC„, 

i n ""(Cife, *<*))- n v"{c 2 (t u t ai ))\ < ni /'(c) B -v"(c)AC( S ). 

2£cr((T) iGcr(cr) 

Thus, with r = f3^v"{c) we easily deduce from (4.9) that 

A#(s,i) < 4/3 2 ^"(c)( S -t) 2 ^n(2n!)- 1 [2r(s-t)] 2 ("- 1 )AC(s) 

n>l 

(4.21) < 2l3 2 iy'"(c)(s-t) 2 e 2r ^AC(s). 

Next, note that ip'(-) is a polynomial of non-negative coefficients, hence non-decreasing on M + . Consequently, 
for any (Ri, d) e -4(c), r, p, c) and all (u, u) € T, 

A(^(C)i2)(«,«) < ^(c)[J?i(»,!i)AC(b,u) + C 2 (v,u)AR(v,u)} < ip'(c)[pAC(v, u) + cAR(v, u)]e~ s{v - u) . 

As \e~ x — e~ v \ < \x — y\ for all x, y > 0, we thus get that 

AA(M) < (3 2 e- (s - t)/2 f f A{^{C)R){v,u)dudv 

Jt Jo 

/S PV 
/ e'^'^dudv 

(4.22) < c 2 l3 2 { S -t)e- (s - t ^ 2 [AC{s) + AR(s)}, 

for c 2 = ip'(c)(p+ c)J _1 . Since R(s,t) = H(s,t)A(s,t) we now obtain from (4.11), (4.21) and (4.22) that 
AR(s,t) < A 1 (s,t)AH(s,t)+H 2 (s,t)AA(s,t) 

< 2/?V"(c)(s - t )2 c -[5-2r]( s -i) AC ( s ) + ClC2/3 2 (s _ t)e-[i- 2r ]('-*)[AC7(«) + Ai2(a)] 
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Further, as 1/2 — 2r > 1/6, the latter bound leads to (4.18) for the finite universal constant 

M R := sup (2v'"(c)9 2 e- d ^ + Cl c 2 9e- 9 ^) . 
e>o 

• The Lipschitz bound (4.19) on C. 

Recall that C(s, t) for (s, t) G T is given by (3.1), whereas for all (R, C) G A(S, r, p, c) and < t < v, 

h(v,t)= C{u,t)R(v,u)u"(C(v,u))du < cpKr^v" ' [c)e- s(v - t) , 
Jo 

I 2 (v,t)= f u'(C(v,u))R(t,u)du < cpKr~ 1 iy"(c)e^ s{v - t) . 
Jo 

Further, with c 3 = p(i>"{c) +cv'"(c)), c 4 = cv"(c), c 5 = <5 _1 max(c 3 , 2c 4 ) and c 6 = <5 _1 i/'(c) max(p, 2c), it is 
not hard to check that for (i? i; d) G A(S, r, p, c) and v > t > 0, 

Ah(v,t) < c 3 AC(v) [ e- 5{v -^du + Ci AR{v) f e -*|«-*l dti < C5 [AR(v) + AC{v)] 
Jo Jo 

AI 2 (v,t) < pv"{c)AC{v) f e- 5( - t -^du + cv"(c)AR(t) [ e'^'^du < c 6 [AR(v) + AC(v)} 

Jo Jo 

Thus, with c 7 = cpK/'S and /?V (c)^ 1 = r < 1/6, we have from (4.22) that for s > t > 0, 

AC(s,t) < AA(s,t)+c 7 [ AA(s,v)dv + /3 2 f e-^ s - v \Ah(v,t) + AI 2 (v,t)]dv 

Jt Jt 

< p 2 M c [AC(s)+AR(s)], 

where 

/>oo 

M c := c 2 sup 0e~ e/2 + c 7 c 2 / ue~ u/2 du + 2(c 5 + c 6 ) , 
e>o Jo 

and consequently, (4.19) holds. 

Suppose 5,r,p,c are such that f is a contraction on A(8, r, p, c), hence also on its non-empty subset 
S(S, r, p, c). Starting at some So = (Ro, Co) G S(S, r, p, c) consider the sequence Sk = ^{Sk-i), k = 1, 2, . . ., in 
S(S, r, p, c). Since ^ is a contraction, clearly {Sk} is a Cauchy sequence for the uniform norm || • || of (4.14). 
Hence, Sk — > Soc in the Banach space (C(T) x C(R + x K + ), j| • ||). Note that A(S, r, p, c) is a closed subset of 
this Banach space, so Soc G .4((5, r, /?, c). Further, fixing r > 0, with |(£,y)| := |a;| + \y\, since Sk € S{8 1 r 1 p 1 c) 
we have that 

lim sup \S 00 (t + T,t)-S QC (t l + r,t')\ < 2||5 0O -5 fc ||+ lim sup \S k (t + r, i) - S k (t' + r, i')l = 2||Soo - Sfc|| • 

T ^°°t,t'>T T ^°°t,t'>T 

Taking fc — > oo we deduce that {Soo(i + r, t)} is a Cauchy function from R + to [0, p] x [0,c], hence 5oo(t + r,t) 
converges as t — > oo. With this applying for each r > 0, we see that 5^, G S(S,r, p,c) and further that S^ 
is the unique fixed point of the contraction ^ on the metric space (S(S, r, p, c), || • ||). By our construction of 
<J>, it follows that the fixed point S^ = (i?, C) of ^ satisfies (1.14)-(1.15), from which we conclude that the 
unique solution of the latter equations is in S(6,r, p, c). As noted before, this yields the existence of R{dt(T~) 
of (4.4) and Cf dt (r) of (4.5), such that (Rmt, Cfdt) € B(S, r, p, c). 
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In the course of proving Proposition 4.1 we found that on S(6,r,p,c) the mapping induces a mapping 
*fdt : (-Rfdt, Cfdt) -» (-Rfdt, Cfdt) such that 

. 2n 

i? fdt (r) = A»]T/? 2 " ]T / J] ^(C fdt (^ - ^ W )) J] <®i - 

™>o CT eNC„ o - 01 -"- e2 "- T »ecr(<T) J=l 

Cfdt(r) = A(r)+/3 2 / A(t - v)h{v)dv + (3 2 [ A{t - v)T 2 (v)dv , 
Jo Jo 

where A(r) = e~^ T for /U = l/2 + /3 2 / of (4.17), while Ti(v) and J 2 (v) of (4.12)-(4.13) arc the two integrals on 
the right-hand-side of (4.16). In particular, i?f dt and C {dt are differentiable on R + , and by [12] we have that 
for r > 0, 

(4.23) R' Mt (r) = -fiR (dt (r) + (3 2 f R {d t(r - 6)^(6)*" {C at (6))d0 , 

Jo 

with .Rfdt(O) = 1, while 

(4.24) C' fdt (T) = -pC m (r) + f3 2 h{r) + /3 2 / 2 (r) , 

with Cfdt(0) = 1. Since the solution (R,C) of (1.14)-(1.15) is a fixed point of the corresponding pair 
(i?fdt,Cf dt ) is a fixed point of * fdt , which by (4.23)-(4.24) satisfies the FDT equations (4.15)-(4.17) with 
1 = 0. 

Recall that B{8, r, p, c) consists of pairs [R, C) of functions that are uniformly bounded, exponentially 
decaying to zero at infinity, and of fixed values at zero. It is not hard to verify that when 1 = 0, any 
solution of (4.15)-(4.17) in this space, is a fixed point of \&f d t- Further, \I/fdt is a contraction on B(5,r,p, c), 
equipped with the suprcmum norm, since following the very same arguments we used in proving the Lipschitz 
estimates (4.18)-(4.19) of Proposition 4.2, we find that Afi fdt (r) < M R f3 2 [AR idt (T)+AC [dt {T)} and AC fdt (r) < 
M c f3 2 [Ai?fdt (t) + ACfdt(r)] for all t < oo, now with A/(t) = sup 0<r \f(0)\. This of course proves the 
uniqueness of the solution of (4.15)-(4.17) in B(5, r, p, c), in case I = 0, as claimed. □ 



5. Study of the FDT equations 

In this section, we complete the proof of Theorem 1.3 by relating the solutions of the equations (4.15)-(4.17) 
with the solution of (1.20). 

Specifically, we first prove Proposition 1.4 about existence, uniqueness, limiting value and exponential 
convergence of the solution of (1.20). Fixing (3 < (3 C , we then know that the unique solution of (1.20) for 
b = 1/2 and 4>{x) = b+2(i 2 v' '(x), exists, is twice continuously differentiable, positive, and of negative derivative, 
such that both converge exponentially fast to zero when s — > oo, with positive exponent ep = | — 2/3 2 z/'(0). 
Further, for [3 < j3\ < (5 C sufficiently small, by Proposition 4.2, the pair of functions (R(t + t,t), C(t + t,t)) 
converges for t — > oo to the unique solution (R{ d t{ T ),Cfdt{ T )) of (4.15)-(4.17), with 1 = 0. Hence, our next 
proposition (whose proof is provided at the end of the section), completes the proof of the theorem by showing 
that for any /3, the solution (Cf d t,-Rfdt) of (4.15)-(4.17) can be expressed in terms of the solution of (1.20) for 
6=1/2 and 4>(x) = 7 + 2(5 2 v' '(x), provided the constant 7 is chosen accordingly. That is, with I 7 matching the 
constant I of (4.17) while (/>(■) satisfies the condition (1.21) of Proposition 1.4. Indeed, if (3 < f3 c , then for 7 = b 
we have that Doo = 0, I 7 = and further, the resulting solution (R, C) has the exponential decay property of 
Proposition 4.2 provided (3 is small enough (for the condition b + 2(3 2 v'(l) > 2^/b2(3 2 v" (\) of Proposition 1.4 
to hold). 

Proposition 5.1. Suppose D{s) is a positive, twice continuously differentiable and decreasing solution of 
the equation (1.20) for b > 0, > and </>(•) which satisfy the conditions of Proposition 1.4. Then, 
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R(s) = -b^D'is) > 0, C(s) = C(-s) = D(s) > and fx = (f)(1) > are such that R(0) = C(0) = 1 and for 
any 7 e R, 

(5.1) R\t) = -nR{ T )+bf R(T-6)R(6)<j>'{C(6))d6, 

Jo 

r-OO r-OC 

(5.2) C'(t) = -nC{r)+bJ C{t - 9)R(6)<f>' (C{6))d6 + b J [<f>(C{0)) - i\R(9 - r)d6 + I 7 

(5.3) /x = 6 + ^ 7 (C(6')) J R(6i)d6' + I 7 , 

Jo 

for V> 7 (x) = x<fi'(x) + 4>{x) — 7 and I 7 = 7 — 6 + .Doo(0(Ax>) — 7) ( r ^call also that b = </>(D 00 )(l — -Doo) 

Proof of Proposition 1.4. Let 7 = [0, 1 — Doo] and M oa (M. + ,I) denote the space of continuous /-valued 
functions on M + which take the value 1 — at time zero. Set <p(x) = 4>{x + -Doc) — 0(1? 00) and for any 
E e Af 00 (R+,7), let k(s) = bip'(E(s)) and 



(5.4) H s (E) = l + J2 e / n ^Mo-^n**- 

n>l cr£JVC n J0<ti---<t 2 n<s , ecr ( CT ) »=l 

Note that the continuously differentiable s 1— > H S {E) is the solution of 

(5.5) = 6 f - v))H s _ v (E)H v (E)dv, H (E) = 1 . 

as Jo 

For /Li = 0(1) let 

$(E)(s) = (l-Doo-6 / e-» v H v (E)dv) V0. 
Jo 

Since 6<p' is non- negative on the interval /, it follows that H S (E) is non- negative and consequently s 1— » $(_E)(s) 
is non-increasing and also belongs to M 00 (R + ,7). Further, since ip' is by assumption Lipschitz continuous on 
the compact set I it follows from (5.4) that for each T < 00 there exists Ki(T) < 00 such that for all 
E, E e and u<T, 

\H U {E) - H U {E)\ < ki(T) f \E(v) - E(v)\dv 

Jo 

With w 1 — > y V Lipschitz continuous, this implies that for some finite, non-decreasing K2(T) > 1, all E, E £ 
M 00 (R+,7) and s < T, 

- $(£)(s)| < «2(T) / |£(w)-£( V )|aV 
Jo 

Thus, $(•) is a contraction on M 00 (M + ,7) equipped with the weighted Li-norm := J °° \E(s)\w(s)ds 

for weight function w(s) — exp(— 2 J Q K2(u)du) > 0. Consequently, we deduce that &(E) = E has a unique 
solution in M 00 (R + ,J), denoted hereafter by E*(s). Let cr* > be the first s > where $(.E*)(s) = 0. Note 
that the function s $(_E*)(s) is continuously differentiable on the interval [0, cr*), in which case also 

E.(s) = 1 - £»oo - b f e-» v H v (E*)dv , 
Jo 

is twice differentiable, with E^s) = —be~ lis H a (E^). Hence, applying (5.5) for H S (E*), we have that for 
s E [0,(J*), 

E';{s) = -fiE'M - f <p'(E.(s - v))El(s - v)K{v)dv. 
Jo 



LIMITING DYNAMICS AT HIGH TEMPERATURE 



19 



Integrating this equation we find that 

K{s) = K(0) - f\fi + <p(E*(s - v)) - tp(E*(0))]E:(v)dv . 
Jo 

Since E' t (0) = -b, E*(0) = 1 - D x and \i - ^(£*(0)) = 0(Ax>). it follows that on [0, cr*) 

(5.6) K(s) = -b-l <j>{E*{s-v)+D 00 )E'«(v)dv, E,(0) = 1 - £)«, . 

Jo 

We next show that cr, = oo, which in view of (5.6) yields that D(s) = E*(s) + D x is a twice differentiablc, 
strictly decreasing [D^, l]-valued solution of (1.20). To this end, note that H S (E^,) > 1 since E 1 * £ M 00 (WL+, I) 
(see (5.4)). Consequently, for all v < cr* both E^(v) < —be~^ v < and 

-4>(E.(s -v) + D^E'M > . 

Thus, from (5.6) we have that for s < cr*, 

(5.7) -be-» s > Ei(s) > -b- 0(£) oo )(^(s) - E.(0)) = -^(D^E^s) , 

since & = ^>(£) 00 )(1 — £>oo) = </ ) (^oo)^'*(0) by the definition of D m . If cr* < oo, then as s f cr* we have by the 
continuity of £?*(•) that £*(s) — > £*(cr*) — <I>(i?*)(cr*) = 0, in contradiction with (5.7). 

The uniqueness of the solution to (1.20) follows from the uniqueness of the preceding £?*(•). Indeed, if a 
continuously differentiablc [D^, l]-valued D(s) solves (1.20), then E(s) = D(s) — £ M 00 (M + , I) and since 
b = (j>(Doo)E(Q) we have that 

(5.8) E'{s) = -(j>{D 00 )E{s)- f ip(E(s-v))E'(v)dv, E(0) = l-D oo . 

Jo 

Hence, E(s) is twice differentiablc, with E'(0) = —b, and differentiating (5.8) we get that 

(5.9) E"(s) = -nE'{s) - f <p'{E{s - v))E'{s - v)E'(v)dv 

Jo 

(using the fact that (fi(E(0)) + (f>(D x ) = fj,). Thus, H S (E) = —b~ 1 e~ tls E'(s) solves (5.5). This means that for 
all s > 0, 



(5.10) E(s) 



= l-D OQ -b f e-» v H v (E)dv . 
Jo 



Since E(s) > 0, so is the right-hand-side of (5.10), that is E = &(E). Since $(•) has a unique fixed point in 
M 00 (R + ,7), the solution of (1.20) must also be unique. 

Next, the monotone and bounded function D converges to some x £ [D^,!]. Fix M — M e such that 
D(v) < x + e for all v > M. Then, for s > 2M we have that 

\ (j){D(v))D'(s- v)dv = / (j>(D(v))D'(s-v)dv+ / 4>{D(v))D' (s - v)dv 
Jo Jo Jm 

> -<j)(l)(D(s - M) - D(s)) - <j>(x + e)(l - D(s - M)) 

> - (x + e - x) - (f){x + e){l - x) . 

Hence, for all s > 2M, 

D'(s) <-b + 4>(l)e + (f>(x + e)(l - x) . 

Since D(s) is bounded below, it follows that (j>(l)e + <p(x + e)(l — x) > b. Taking fjOwe see that x = 
(since x > and is the largest y £ [0, 1] for which <fi(y)(l — y) > b). 
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Finally, recall that the function E(s) = D(s) — -Doc is strictly positive, monotone decreasing with E'(s) = 
-be~^ s H s {E) for H S (E) of (5.4). Since there k(s) = btp'{E{s)) < b<j>'(l) (by the assumed convexity of </>(•)), 
the same argument as in (4.10) yields that for all s > 0, 

/2 
e ^/wwsx^ A _ x 2 dx < Cie v^'W s . 

Consequently, if 0(1) > 2- N /b0'(l), then both E'(s) and D(s) converge exponentially to zero for s — > oo. □ 
Remark. We note in passing that <p(x) > <p'(Q)x for all a; G [0, 1] when is convex, in which case by the 
monotonicity of E(s), 

J(s) := ^(s)- 1 / - v))dv > (a;-V)(-B(*))[^(0) - £(*)] > ^(0)[1 - Ax, - £(*)] • 

Jo 



/o 

Further, recall that from (5.8), 



<,]0 " E (s) = %(s) = J(s)-4>(D 00 ). 



ds E 

hence in case <j){D oo) =^'(£) 0O )(l-£) oo ),wefind that £(s) > S(0) cxp(-< ? !)'(L> 00 ) / Q s E(u)du), which thus does 
not converge to zero exponentially fast as s — > oo. 

Proof of Proposition 5.1. Since D(0) = 1 and D(9) J. Doo, for i? = —b~ 1 D' and any 7 € R we have that 

roc 

76 / R{6)d6 = 7(1 - Doo) = I 7 - Io , 
Jo 

hence we may and shall assume hereafter that 7 = 0. Further, it is easy to check that for ipo{ x ) = \x<j>{x)]' 

/>oc 

/ MD(6))D'(6)d0 = £>«,#£>«,) - D(O)0(£>(O)) , 
Jo 

so the choices of \i = 4>{1) = </>(D(0)) and I = — b + £ , oc (/'(£ ) oc) guarantee that C — D and R = —b~ 1 D' 
satisfy (5.3). Recall that while proving Proposition 1.4 we have seen that D'(0) = —b and for all r > 



D"(t) = -hD'{t) - f D'(t - 6)D'(8)<j>'(D(0))d0 
Jo 



(c.f. (5.9)), from which it immediately follows that C = D and R = —b~ 1 D l satisfy (5.1). 

Turning to the remaining task of verifying that (5.2) holds, note that since D(-) is bounded and converges 

to Doo, 



Similarly, 



/oo poo 
D(9 - T)D'(6)cj>'{D(6))d6 + J cj>(D(6))D' {6 - r)d8 = <KAx,)Ax, - 0(£>(t))D(O) . 

f D{T-6)D'{6)<j>'{D{6))d6 - f 4>(D(0))D'(t - 6)d6 = <P(D{t))D{0) - 4>(D(0))D(t) . 
Jo Jo 



Hence, using the symmetry C(s) = C(—s), upon substituting in (5.2) C = D, R = —b 1 D' and the preceding 
two identities, it is not hard to verify that (5.2) holds for 7 = if 



D'(t) = -0(Doo)Doo + Io _ / <t>( D ( e )) D ' ( T - e ) d0 > 

Jo 



which in view of (1.20) is merely the statement that b = </>(£> oo)!?^ — Io. Our choice of Io guarantees that 
the latter identity applies, thus completing the proof of the proposition. □ 
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6. About the FDT regime at all temperatures 

We return to the equations (1.15)— (1.16) and explain why (1.20) is the natural candidate for describing the 
function Cfdt(-) for all (3. To this end, set for s > t, 

G{s,t) = R(s,t) - 2d t C{s,t) 

and 

I(s, t) = <3 2 [ [C(t, u)G(s, u)v"{C(s, u)) + u'(C(s, u))G(t, u)]du - 2/?V(C(s, 0))C(t, 0). 
Jo 

Since ip(x) = [xv'(x)]' and C(s, s) = 1, it follows that 

I(s,s) = (i 2 f i/j(C(s,u))G(s,u)du-2p 2 u'(C(s,0))C(s,0) 
Jo 

f 

Jo 

for p of (1.16). Further, by similar reasoning, 



((). !) = [3 2 I iP(C(s,u))R(s,u)du-2f3 2 v'(l) = p-^-20 2 v'{\) 



2/?V(C(s,i)) +/(s,i) = 2(3 2 / d u {C(t, u y(C(s,u))]du + 2l3 2 i>'(C(s,0))C(t,0) + I(s,t) 
= 



t 

2 



/ [C(t, u)R(s, u)v"{C(s, u)) + u'(C(s, u))R(t, u)\du . 
Jo 



Thus, in these notations (1.14)-(1.16) are equivalent to having for s > t, 

(6.2) d s R(s,t) = -[p + I(s,s)]R(s,t) + (3 2 £ R{s,u)R{u,t)v" '(C(s,u))du , 

(6.3) d s C(s,t) = -[ P + I(s lS )}C(s,t)+f3 2 J C(u,t)R(s,uy'(C(s,u))du + 2/3 2 v'(C(s,t))+I(s,t), 

with R(t, t) = C(t, t) = 1 and p = 2" 1 + 2/?V(l). 

We note in passing that G(s, t) is such that for s > t, 

(6.4) d,G(a,t) = -2d t I{s,t)-[p + I(s,s)]G{s,t) + 2{3 2 v"(C{s,t))G(s,t) 

+ f3 2 J G{u,t)[G{s,u)u"{C{s,u) + 2d u C{s,u)u"{C{s,u))}du. 

The physics prediction is that G(t + r, t) — > as t — > oo while r is fixed, for any finite (3 (this is the famous 
FDT relation). As a result, comparing I(s,t) and I(t,t), we further expect that I(t + r, t) — * I when t — > oo 
while r is fixed. We next show that the latter ansatz results with the existence of an FDT solution (Rfdt, Cfdt) 
such that Rf At = -2C' m and C fdt solves (1.20) for b = 1/2 and <p(x) = 7 + 2/?V(x), where 7 = 1+ 1/2. 

Proposition 6.1. Assume that given the continuous function I(s,t) there exists a continuously differ entiable 
solution (Ri(s,t),Ci(s,t)), s > t > 0, of (6.2)-(6.3) with p = 2~ x + 20 2 v'(l), the initial conditions R(t,t) — 
C(t,t) = 1 and uniformly bounded C(s,t). Further, suppose that for any T < 00 



(6.5) lim sup |7(t + r,i)-I| =0, 

*^°°tG[0,T] 

where the constant I is such that 

(6.6) sup {{l+\ + 2[3 2 V '{x))(l-x)}> 1 - 
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Then, Ci(t + r, t) — > Cfdtfr) as f — > oo, uniformly in t S [0, T], where Cfdt(T) is i/ie unique solution of (1.20) 
forb=l/2 and(j)(x) = j + 2(3 2 i/ (x), with 7 = 1+1/2. Further, (i + r, i) -» ifot (r) = -2Cf dt (r) as t -> 00,. 
uniformly in t G [0, T] . 

Remark. In Theorem 1.3 we circumvent the difficulty of showing the ansatz that J(t + r, t) converges to a 
(/3-dcpendent) constant at any value of f3, by verifying that at very high temperature, i.e. sufficiently small 
[3, the exponential decay to zero in s — t of (R, C) results with the convergence to zero of /. Proposition 6.1 
shows that this is the only obstacle to extending our results about the FDT regime to all temperatures. 

Before proving the proposition, consider its consistency with our choice of j((3) for (3 > (3 C based on 
q([3) of (1-24). Specifically, considering in (6.1) the contribution to the integral from u e [s — M,s] with 
M < 00 arbitrarily large, we expect the FDT solution (i?fdt,Cfdt) to contribute — 2(3 2 q{[3)i>' {q([3)) to the 
limiting constant I = j(/3) — 1/2 (for which (6.6) holds by our choice of j((3))- We further expect I to be the 
sum of this FDT contribution and a non-negative contribution from the aging regime (i.e. the integral over 
u £ [0, s — M\). Given the relation between p and I 7 in Proposition 5.1, we deduce that I 7 is exactly the 
contribution of the aging regime to I. In conclusion, we should have I 7 = for (3 < f3 c , I 7 > at (3 = [3 C and 
I 7 > for (3 > (3 C . This is indeed the case, for when (3 < (3 C we have 7 = b = 1/2 and = leading to 
I 7 = while for (3 > (3 C , since v'{x) < xv" (x) for x > with strict inequality when x > 0, our choice of 7(/3) 
leads to 

2/ 7 = 4/3 2 (l - q((3W'W)) - Aq(0))] - 1 > 4/3 2 (l - q{(3)fv" (q([3)) - 1 > , 
with a strict inequality whenever q{J3) > (in particular, for all (3 > (3 C ). 

Proof of Proposition 6.1. Fixing hereafter the values of (3 and p = l/2 + 2/3 2 j/(l) and the initial conditions 
R(t, t) = C(t, t) = 1, a re-run of the argument at the end of the proof of Proposition 1.1 shows that per given 
continuous function I(s,t), the system of equations (6.2)-(6.3) admits at most one bounded solution, denoted 
(Ri(s, t), Ci(s, t)) on any compact interval < t < s < T. Consequently, it has at most one continuous 
solution (Rj, Ci) for all s > t > 0. 

In particular, in case I(s,t) — I is a constant that satisfies (6.6), we know that (1-21) holds for b = 1/2, 
7 = 1 + 1/2 and 4>{x) = 7 + 2(3 2 v'{x). With D(t) denoting the unique solution of (1.20) for these parameters, 
we claim that -Rj(s, t) = —b~ 1 D'(s — t) and Gj(s, t) = D(s — t) is then a solution of (6.2)-(6.3) (hence its unique 
solution). Indeed, recall that D(0) = 1, D'(0) = —b so the given initial conditions for (i?j,Gj) hold. Further, 

these settings result with p := p + I = 0(1), so upon taking u ~ t + 9 and s = t + r we see that the equation 
(6.2), which is a direct consequence of (5.1) of Proposition 5.1, clearly holds. By the same transformations we 
see that our proposed solution satisfies (6.3) provided that 

(6.7) D'(t) = -c/){1)D(t) - [ D(t -9)D'(9)(t)'(D(9))d9 + 4>(D(t)) -b. 

Jo 

Noting that by integration by parts, 

/ D(T-9)D'(9) ( f>'(D(9))d9 = -D(T-9) ( p(D(9))\ T Q + [ D'(t - 6)<j>{D(9))dB , 
Jo Jo 

and having D(0) = 1, we see that (6.7) is equivalent to (1.20) hence holds as well. 

Fixing T < 00, for each £ positive let :— {(s, t) : £ < t < s < t + T} C R+ x R + . We proceed by showing 
that per fixed £ > and T < 00 the mapping I (i?/, Ci) is Lipschitz with respect to the supremum norm 
over with a Lipschitz constant that is independent of £. More precisely, given two bounded functions I(s, t) 
and T(s,t) to which correspond continuously differentiable solutions (Ri, Ci) and (i?/, Cj) of (6.2)-(6.3), such 
that Ci and C/ are uniformly bounded, we let A^I :— sup{|/(s, t) — I(s, t)\ : (s, t) e T^}, and show that 



(6.8) 
(6.9) 



A e R := sup{\R I (s,t)-R T (s,t)\:{s,t)GT i }<K A i I, 
A^C := sup{\C I (s,t)-C I (s,t)\:(s,t)eT i }<K A i I, 
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where kq is a finite constant that depends only on T, /?, v' , p and the uniform bound on 7, 7, Cj and 
Cj for (s,t) e r^. Considering the bounded I = I, we have seen already that Rj(t + r,t) = Rfdt(T) and 
Cj(t + r,t) = Cfdt(T) with Cj bounded. Our assumption (6.5) then amounts to A^7 — > as £ — > oo. In 
particular, this implies that 7 is also uniformly bounded on T^, so from (6.8) and (6.9) we get that A^Ti. — > 
and A^C — ► 0, which are easily seen to match our desired conclusion. 

We thus complete the proof by verifying the bounds of (6.8) and (6.9), by an argument similar to the 
one we used for proving uniqueness of the system of equations (1.14)-(1.16). To this end, note first that 
Ri = A/T/cj for Aj(s,t) = exp(— J t (p + I{u,u))du) and Hc(s,t) of (2.2). Thus the assumed uniform bounds 
for 7 and C/ on imply that Ri is also uniformly bounded on by a constant that depends only on 
T, (3, p, v" and the corresponding uniform bounds for 7, Cj. Of course, the same applies for Rj. Next 
let AI(v,u) = \I(v,u) - I(v,u)\, AR(v,u) = |i?j(u,u) - Rj(v,u)\, AC(v,u) = \d(v,u) - Cj(v,u)\ and for 

(s,t) e r 5 set 



h(s,t) 



J [AR(s,u) + AC{s,u)]dv 



Then, similarly to the derivation of (2.7) and (2.8), upon considering the difference between the integrated 
form of (6.2)-(6.3) for our solutions (Cj, Ri) and (C/, Rj), we find that for any (s, t) G 

(6.10) AR(s,t) < kx[[ AR(v,t)dv+ [ h{v,t)dv + f AI(v,v)dv], 

Jt Jt Jt 

(6.11) AC(s,t) < naif AC ( v >t) dv + [ h(v,t)dv+ ( AI(v,t)dv+ [ AI(v,v)dv], 

Jt Jt Jt Jt 

where the positive K\ < oo depends only on T, (3, p, u'(-) and the maximum of 7, |7?/|, |C/|, |7|, \Rj\ and \Cj\ 
on r^. Replacing t by u, summing these two inequalities and then integrating the result over u e [t, s] yields 
that 

(6.12) h(s,t)<K 2 [ h(v,t)dt+ AI(v,v)dv+ / AI(v, u)dudv] , 

Jt Jt Jt Jt 

for any (s, t) € Tj, with a finite, positive constant k 2 (of the same type of dependence as K\). Since h(t, t) = 
for all t, we get from (6.12) by Gronwall's lemma that h(s,t) < k^A^I for some finite, positive K3 (of same 
dependence type as K\) and all (s,t) £ T^. Recall that AR(t,t) — AC(t,t) = due to the given initial 
conditions, so upon plugging into (6.10) and (6.11) our uniform bound on h(s,t), we complete the proof of 
(6.8)-(6.9) by yet another application of Gronwall's lemma. □ 
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